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Abstract. This paper develops some interior penalty hp-discontinuous Galerkin 
(hp-DG) methods for the Helmholtz equation in two and three dimensions. 
The proposed hp-DG methods are defined using a sesquilinear form which is 
not only mesh-dependent but also degree-dependent. In addition, the sesquilin- 
ear form contains penalty terms which not only penalize the jumps of the 
function values across the element edges but also the jumps of the first order 
tangential derivatives as well as jumps of all normal derivatives up to order 
p. Furthermore, to ensure the stability, the penalty parameters are taken as 
complex numbers with positive imaginary parts. It is proved that the proposed 
hp-discontinuous Galerkin methods are absolutely stable (hence, well-posed). 
For each fixed wave number k, sub-optimal order error estimates in the broken 
.ff^-norm and the L 2 -norm are derived without any mesh constraint. The er- 
ror estimates and the stability estimates are improved to optimal order under 
the mesh condition fc 3 h 2 p _1 < Co by utilizing these stability and error esti- 
mates and using a stability-error iterative procedure To overcome the difficulty 
caused by strong indcfinitcncss of the Helmholtz problems in the stability anal- 
ysis for numerical solutions, our main ideas for stability analysis are to make 
use of a local version of the Rellich identity (for the Laplacian) and to mimic 
the stability analysis for the PDE solutions given in 1 191 1201 33 , which enable 
us to derive stability estimates and error bounds with explicit dependence on 
the mesh size h, the polynomial degree p, the wave number fc, as well as all 
the penalty parameters for the numerical solutions. 



1. Introduction 

This is the second installment in a series (cf. [27]) which devotes to developing 
and analyzing novel interior penalty discontinuous Galerkin (IPDG) methods for 
the following Helmholtz problem with large wave number: 

(1.1) -Au-k 2 u = f inQ:=fii\D, 

(1.2) h iku = g on T R := dtt u 

(1.3) u = onT D :=dD, 
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where k E H, called wave number, is a (large) positive number, D C fii C R d , d — 
2, 3, D is known as a scatterer and is assumed to be a bounded Lipschitz domain, 
f2i, which is assumed to be a polygonal/polyhedral domain and often taken as a 
(i-rectangle in applications, defines the size of the computational domain. Note that 
dVl = Tr U Tjj. nfi denotes the unit outward normal to O. i := \f— 1 denotes the 
imaginary unit. Condition (| 1 . 2|) with g — is known as the first order absorbing 
boundary condition (cf. |23j). which is used to minimize the reflection at the 
boundary Tr and to limit the computation of the original scattering problem just 
on the finite domain O. Boundary condition (|f .3[) implies that the scatterer is 
sound-soft. We note that the case D = also arises in applications either as 
a consequence of frequency domain treatment of waves or when time-harmonic 
solutions of the scalar wave equation are sought (cf . [22 ) . 

In [37] we proposed and analyzed some IPDG methods for problem (jl.l|) - (|1.3|) 
using piecewise linear polynomial trial and test functions. It was proved that the 
proposed methods are unconditionally (with respect to mesh size h) stable and 
well-posed for all wave number k > 0. Optimal order error estimates were estab- 
lished showing explicit dependence of the error bounds on h, k and all penalty 
parameters. However, due to the existence of a pollution term, the (broken) H 1 - 
norm error bound deteriorates as the wave number k increases under the practical 
"rule of thumb" mesh constraint kh < 1. To improve the accuracy and efficiency 
of those IPDG methods, it is necessary to use (piecewise) high order polynomial 
trial and test functions partly because of the rigidity and low approximability of 
linear functions and partly because of the very oscillatory nature of high frequency 
waves. However, simply replacing the linear element by high order elements in the 
IPDG methods of [37] does not reduce the pollution very much, in particular, the 
theoretical error bounds do not change much because the analysis of [27] indeed 
strongly depends on the properties of linear functions. 

Motivated by the above challenge and observation, the primary goal of this 
paper is to develop some new /ip-interior penalty discontinuous Galerkin (hp-IPDG) 
methods which retain the advantages of the IPDG methods of [57J but improve their 
accuracy and stability by exploiting the efficiency and flexibility of piecewise high 
order polynomial functions. To the end, our key idea is to construct a sesquilincar 
form (as a discretization of the Laplacian) which is not only mesh-dependent (or 
/i-dependent) but also degree-dependent (or p-dependent) by introducing penalty 
terms which not only penalize the jumps of the function values across the element 
edges but also the jumps of the first order tangential derivatives as well as jumps of 
all normal derivatives up to order p. In addition, as in 27J, to ensure the stability, 
all penalty parameters are taken as complex numbers with positive imaginary parts. 
Since the Helmholtz equation with large wave number is non-Hermitian and strongly 
indefinite, as expected, stability estimates (or a priori estimates) for numerical 
solutions under practical mesh constraints is a difficult task to carry out regardless 
which discretization method is used. To overcome the difficulty, as in [27], the 
cruxes of our analysis are to establish and to make use of a local version of the 
Rellich identity (for the Laplacian) and to mimic the stability analysis for the PDE 
solutions given in [T9l [20 , 33J. The key idea here is to use the special test function 
Vuh -(x — xn) (defined element-wise) with Uh denoting the /ip-IPDG solution, such 
a test function is valid for any DG method. We remark that the same technique 
was successfully employed by Shen and Wang in [38] to establish the stability and 
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error analysis for the spectral Galerkin approximation of the Hclmholtz problem. 
We also note that although the similar techniques to those in [2TJ, [38] are utilized 
in this paper to carry out the stability analysis, the analysis of this paper is more 
involved because the special sesquilinear form of this paper, which contains jumps 
of high order normal derivatives, is a lot more complicate to deal with, even they 
are similar conceptually. 

Since the Helmholtz equation appears, in one way or another, directly or indi- 
rectly, in almost all wave-related problems arisen from many science, engineering, 
and industry applications, solving the Helmholtz equation, in one form or another, 
has always been and remains at the center of wave computation. We refer the 
reader to (P3 El ISl IH3 U31 H3I U51 HS1 1221 1231 12SI 1321 stn<i tlae refer- 
ences therein) for some recent developments on numerical methods, in particular, 
Galerkin type methods, for the Helmholtz equation. We also refer the reader to [27] 
for a brief review about some theoretical issues for finite element approximations 
(and other types of Galerkin approximations) of the Helmholtz equation. 

The /ip-finite element method (/ip-FEM) is a modern version of the finite ele- 
ment method, capable of achieving exceptionally fast (exponential) convergence. It 
combines the flexibility of the standard finite element method and the high order 
accuracy of the spectral method. Consequently, the /ip-FEM can often attain more 
accurate results than the standard finite element method does while using less CPU 
time and resources. The /ip-FEM has undergone intensive developments both on 
theory and implementation in the past twenty five years. We refer the reader to 
the survey paper [7] and two recent monographs (HI 142] for a detailed exposition 
on the basic theory and advanced topics of the hp-FEM. 

Discontinuous Galerkin (DG) methods was first proposed in 1970s, they were not 
popular then because they produce larger algebraic systems than standard finite 
element methods do. However, due to the emergence of high performance com- 
puters and fast solvers since early 1990s, especially, massively parallel computers 
and parallel solvers such as multilevel and domain decomposition methods, which 
together with advantages of DG methods has quickly attracted renewed interests 
in DG methods. They have been heavily developed and tested in the past fifteen 
years, we refer the reader to [4] and the references therein for a review of recent de- 
velopments. As is well known now, DG methods have several advantages over other 
types of numerical methods. For example, the trial and test spaces are easy to con- 
struct, they can naturally handle inhomogeneous boundary conditions and curved 
boundaries; they also allow the use of highly nonuniform and unstructured meshes, 
and have built-in parallelism which permits coarse-grain parallelism. In addition, 
the fact that the mass matrices are block diagonal is an attractive feature in the con- 
text of time-dependent problems, especially if explicit time discretizations are used. 
Moreover, as proved in [27], DG methods are also effective and have advantages over 
finite element methods for the strongly indefinite Helmholtz equation, which has not 
been well understood before. We refer the reader to [3j fl ^ flTl [211 [26l [40], [45] 
and the references therein for a detailed account on DG methods for coercive elliptic 
and parabolic problems, and to [43] [251 EH O [34] and the references therein for 
recent developments on /ip-discontinuous Galerkin (hp-DG) methods. 

The remainder of this paper is organized as follows. In Section [21 we first in- 
troduce notation and gather some preliminaries, and then formulate our /ip-IPDG 
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methods. Both symmetric and non-symmetric methods are constructed and vari- 
ous possible variants are also discussed. Section [3] devotes to the stability analysis 
for the hp-lPDG methods proposed in Section O It is proved that the proposed 
/ip-IPDG methods are stable (hence well-posed) without any mesh constraint. In 
Section [4l using the stability results of Section [3] we prove that for each fixed wave 
number k, sub-optimal order (with respect to h and p) error estimates in the bro- 
ken ii^-norm and the L 2 -norm are derived without any mesh constraint. Finally, 
using the stability estimate of Section [31 the error estimates of Section 0] and a 
stability-error iterative procedure we obtain some much improved (optimal order) 
stability and error estimates for the hp-lPDG solutions under the mesh condition 
k 3 h 2 p^ x < Cq in Section [H where Co is some constant independent of k, h, p, and 
the penalty parameters. 

2. Formulation of Zip-interior penalty discontinuous Galerkin 

METHODS 

2.1. Notation and preliminaries. The space, norm and inner product notation 
used in this paper all are standard, we refer to |H]|15|[9] for their precise definitions. 
On the other hand, we note that all functions in this paper are complex- valued, 
so the familiar terminologies such "symmetric/non-symmetric" and "bilinear" are 
replaced respectively by terms "Hermitian/non-Hermitian" and "sesquilinear" . For 
a complex number a = a r + ia.i (a r and ai are real numbers), a := a r — iai denotes 
the complex conjugate of a. (•, -)q an d (v)s f° r ^ C dQ denote the complex 
L 2 -inner product on L 2 (Q) and L 2 (T,) spaces, respectively. (•,•) := (-,-)o and 
(•, •) := (•, -)qq. We also define 

H^ D (Q) := {u G if 1 (ft); u = on T D ) . 

Throughout the paper, C is used to denote a generic positive constant which is 
independent of k, h, p, and the penalty parameters. We also use the shorthand 
notation A < B and B > A for the inequality A < CB and B > CA. A ~ B is for 
the statement A < B and B < A. 

We now give the definition of star-shaped domains. 

Definition 2.1. Qc R d is said to be a star-shaped domain with respect to xq G Q 
if there exists a nonnegative constant cq such that 

(2.1) (x - Xq) ■ HQ > Cq Vx G DQ. 

Q C R d is said to be strictly star-shaped if cq is positive. 

In this paper, we assume that Qi is a strictly star-shaped domain. Recall that 
ill is often taken as a d-rectangle in practice. We also assume that the scatterer 
D is a star-shaped domain, without loss of the generality, with respect to the same 
point xn x as Qi does. This implies that xq -1 G D C Ox- Under these assumptions, 
there hold following stability estimates for problem (|l.l|) - (jl.3|) . 

Theorem 2.1. Suppose Q,\ C R rf is a strictly star-shaped domain and D C fii is 
a star-shaped domain. Then the solution u to problem (|l.l|) - (|1.3j) satisfies 

(2-2) H| H , (n) < (~ + F"- 1 ) (||/|| L2(0 ) + h\\mv R) ) 
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for j = 0,1 if u e H 3 / 2+£ (fl) for some e > 0. (J2J2|) also holds for j = 2 if 
u G H 2 (fl). Furthermore, there hold 

j-2 

(2-3) \\u\\ H i (Q) < (- + F- 1 ) (||/|U 2 (n) + \\9\\mr R ) + £ ll/ll^ < n) ) 

/or j = 3, 4, • • ■ , q if u G H 2 (fl) n H^ oc {Vl) for some positive integer q > 3. 

Proof. Inequality ([22]) for j = 0,1,2 was proved in [13 [20l [33]. Inequality (f2T3|) 
follows from ([2.2ft and an application of the standard cutoff function technique 
together with an induction argument. We leave the derivation to the interested 
reader. □ 

2.2. Formulation of hp-TPDG methods. To formulate our hp-lPDG methods, 
we need to introduce some notation, most of them were already appeared in [27] , 
Let Tf, be a family of partitions of the domain £1 := f2j \ D parameterized by 
h G (0, ho). For any K G Ti, we define Hk '■— diam(if). Similarly, for each 
edge/face e of K € Th, h e :— diam(e). We impose the following mild restrictions 
on the partition 7^: 

(i) The elements of Th satisfy the minimal angle condition, 

(ii) T is locally quasi-uniform, that is if two elements K and K' are adjacent 
( i.e., me&s(dK n dK') > ), then Hk — hie'- Where meas(e) stands for 
(d — l)-dimensional Lebesgue measure of e. 

For convenience, we assume diam(fi) ~ 1, hence h e ,hx ^ 1. 

For any two elements K, K' e T^, we call e = dK n dK' an interior edge/ face 
of T if meas(e) > 0. Note that e could be portion of a side/face of the element K 
or K' in the case of geometrically nonconforming partition. Also, for any element 
K G T, we call e = dK n d£l a boundary edge/face if meas(e) > 0. Then we define 

Eh '■= set of all interior edges/faces of Th, 

Eh := set of all boundary edges/faces of Ti on Tr, 

E® := set of all boundary edges/faces of Ti on Tp, 

Eh D ■= E^ U Eh = set of all boundary edges/faces of T%, 

Eh D := E^ U E^ = set of all edges/faces of Ti except those on Tr, 

Eh ■= Eh U Eh D — set of all edges/faces of Ti- 

We also define the jump [v] of v on an interior edge/face e = dK n dK' as 

Ml . _ f v \k ~ v \k> , if the global label of K is larger, 
\ v\k' — v\k , if the global label of K' is larger. 

If e G £P , set [u]| e = v\ e . The following convention is adopted in this paper 

{v}\ e :=^(v\ K + v\ K/ ) tfe = dKndK'. 

If e G £f D , set {u}| e = u| e . For every e = dK n dK' G ££, let n e be the unit 
outward normal to edge/face e of the element K if the global label of K is bigger 
and of the element K' if the other way around. For every e G E^ D , let n e = nn 
the unit outward normal to dfl. 
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Let p > 1 be a fixed integer, which will be used to denote the degree of the hp- 
IPDG methods in this paper. For each integer < q < p, we define the "energy" 
space 

E q := \\ H q+1 {K) 7 
KeT h 

and the sesquilinear form a q h (-, •) on E q x E q 

q 

(2.4) a 9 h (u,v) :=b h (u,v) + i(L 1 (u,v) + Y l J j( u , v )) Vw^eB', 



(2.5) b h (u,v):= ^(V U ,V«)k- E ((\ir-}M 

KeT h ee£ in V \ L ° n e J 



d-l 



Pl,eP 



(2.6) Li(u,v) := E /, 

10,eP 



+a ( [u , 



du 



(2.7) Mu,v):= Y, 



dri 

(M,H) e , 



dv 
dv 



dri 



ee£. 



(2.8) J j ( u , v ):=Yl j ,e(j 



ee£. 



dh 



dnl 



3 = 1,2,- 



,q, 



and (J is a real number. jo. e , ■ ■ ■ , 7 9 , e > and (3\. e > are numbers to be specified 
later. {t^}j~} denote an orthogonal coordinate frame on the edge/face e G £h, 
■§^z := Vw • Tg stands for the tangential derivative of u in the direction , and 
denotes the jth order normal derivative of none. 

It is easy to check that (— Au,v) — a q h (u,v) for all u € H q+1 (Vt) and v € 
E q . Hence, a q h (-,-) is a consistent discretization for —A. When a = 1, a q h (-,-) is 
symmetric, that is, a q h {u,v) = a^(w,u). On the other hand, when a ^ 1, a^(-, •) 
is non-symmetric. In particular, er = — 1 would correspond to the non-symmetric 
IPDG method studied in [40] for coercive elliptic problems. In this paper, for the 
ease of presentation, we only consider the case a = 1. The penalty constants in 
(it, v) + Jq(u, v) + ■ ■ ■ + J q (u, w)) are i/3i, e , 170, e, ■ • • , ijq,e, respectively. So they 
are pure imaginary numbers with positive imaginary parts. It turns out that if any 
of them is replaced by a complex number with positive imaginary part, the ideas of 
the paper still apply. Here we set their real parts to be zero because the terms from 
real parts do not help much (and do not cause any problem either) in our analysis. 

Next, we introduce the following semi-norms on the space E q : 



(2.9) 



\l,h 



(E n v < 



Ker h 
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(2.10) 



\l,h,q 



\l,h 



E 



lO.eP 



ll^(e) 



E Pl,eP 
h e 

1=1 



dv 



L2( e ) 



E E 



2j-l 



9 J w 



L 2 (e 



(2.11) 



|l,/>,« 



I l,h,q 



^ D 70,eP 



9u 



L 2 (e 



|i,fc,g are norms on E q if 91? ^ but only semi-norms if 



Clearly, \\-\\ ljh<q and 
<9Z? = 0. 

It is easy to check that the sesquilinear form a*(-, •) satisfies: For any v £ E q 



(2.12) 



(2.13) 



l.h 



Rea^(u,u) = ml 
Ima^(j),«) = Li(v,v) + Jq(v,v) + • • • + J q (v,v). 



Using the sesquilinear form &?(•, •) we now introduce the following weak formu- 
lation for (fn] H [L2] ) : Find u € £? 9 n (fi) n #^(0) such that 



(2.14) 



a? (u, u) - k 2 (u, v) + ik(u, v)r n = (/, v) + (g, v)j 



for any v £ E q n i?p (fi) n il^^O). The above formulation is consistent with the 
boundary value problem (jl.l|l - (|1.2|) because a q h (-, ■) is consistent with —A. 

For any K £ Th, let V P {K) denote the set of all polynomials whose degrees do 
not exceed p. We define our ft,p-IPDG approximation space V£ as 



n v p( r )- 

K<£T h 



Clearly, V,f C E q C L 2 {fl). But V,f ^ ff 1 ^). We are now ready to define our 
hp-IPDG methods based on the weak formulation (|2.14[) : For each < q < p, find 



u\ £ V£ such that 



(2.15) a q h (u q h ,v h ) - k 2 {u q h ,v h ) + ik(u q h , v h )r R = (f,v h ) + (g,v h )r R Vv h £ V%. 



Remark 2.1. (a) When p = q = 1, the above method (|2.15p is exactly the scheme 
proposed in [27]. The L\ term, which penalizes the jumps of the first order tangen- 
tial derivatives, plays an important role for getting a better (theoretical) stability 
estimate in |27j . However, our analysis to be given in the next section suggests that 
the L\ term plays a less pivotal role for high order IPDG methods. 

(b) In fact, (|2.15p defines p+ 1 different IPDG methods for q = 0, 1, • • • ,p. q = 1 
would correspond to using high order elements in the IPDG formulation proposed 
in [27]. 

(c) The idea of penalizing the jumps of normal derivatives (i.e., the J\ term 
above) for second order PDEs was used early by Douglas and Dupont }21j in the 
context of C° finite element methods, by Baker [9] ( with a different weighting, also 
see [26j ) for fourth order PDEs. The idea of using multipenalties Jo, J\,- ■ ■ , J p with 
positive penalty parameters was first used by Arnold in [3] for coercive elliptic and 
parabolic PDEs. The use of L\ term was first introduced in [27] . 
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In the next two sections, we shall study the stability and error analysis for the 
hp-lPDG method (12.15|) . Especially, we are interested in knowing how the stability 
constants and error constants depend on the wave number k (and mesh size h and 
element degree p, of course) and on the penalty parameters, and what are the 
"optimal" relationship between mesh size h and the wave number k. 



3. Stability estimates 



The goal of this section is to derive stability estimates (or a priori estimates) for 
schemes (|2.15|) . To the end, momentarily, we assume that the solution u q h to ()2.15|) 
exists and will revisit the existence and uniqueness issues later at the end of the 
section. We like to note that because its strong indefiniteness, unlike in the case 
of coercive elliptic and parabolic problems (cf. [31 [H [51 [HI HS1 HOI HS]), the well- 
posedness of scheme (|2.15[) is difficult to prove under practical mesh constraints. 

To derive stability estimates for scheme (|2.15[) . our approach is to mimic the 
stability analysis for the Helmholtz problem |TI])-CL2|) S iven in PS EE [33] . The 
key ingredients of our analysis are to use a special test function Vh = a-Vu q h (defined 
element-wise) with a(x) := x — xn t in (|2.15p and to use the Rcllich identity (cf. 
[20] and below) on each element. Due to existence of multiple penalty terms in 
a? h (-, ■), which do not appear in [HI [20l [33] , the analysis to be given below is much 
more delicate and complicate than those of [19l [20l [33] , although they are similar 
conceptually. Since most proofs of this section are in the same lines as those of the 
proofs in Section 4 of [2 7) . we shall omit some details if they are already given in 
[27] . but shall provide them if there are meaningful differences. 

We first cite the following lemma which establishes three integral identities and 
play a crucial role in our analysis. A proof of the lemma can be found in [27] 
Lemma 4.1]. 



Lemma 3.1. 

hold 



Let a(x) :— x — Xn lt v e E 1 , K,K' G 7/j and e € '. Then there 



(3.1) d\\v\\ 2 L 2 (K ) + 2Re(v,a ■ Vv)k = / a-n^W , 

JdK 

(3.2) (d-2)||Vu||i 2m + 2Re(Vu,V(a-Vu))„= / a-n K \S7v\ 2 , 

JdK 



(a ■ n e {Vw} , [Vu]> 



E 

i=i 



dv 1 f dv 1 



dv 



where xq ± denotes the point in the star-shaped domain definition for fli (see Defi- 
nitionWA\). 



Remark 3.1. The identity (f3T2l 
identity for the Laplacian A ( cf. 
for any function » = »/,£ V£ . 



can be viewed as a local version of the Rellich 
321 [20] j. Since V% C E 1 , hence, (ETTj) -rf3~3l) hold 



We also need the following trace and inverse inequalities (cf. [4"T1 14"41 IT^] V 
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Lemma 3.2. For any K £ Th and z e V P (K), 



\\ z \\ L *{dK) ^P h * \\4l'[K) 



Now, taking u/, = u? in (|2.15|) yields 



(3.4) o«(u« u')-* 2 "..'" 2 



«/.:r. + ifc|K" 2 



/i"L 2 (fi) 



Therefore, taking real part and imaginary part of the above equation and using 
(|2.12p and (|2.13[) we get the following lemma. 

Lemma 3.3. Let u q h £ V% solve ([2J5]) . Then 



(3.5) lit 



(3-6) £ 



E ({^}'[<]) e - fc2 KI 



2 

L 2 (n) 



10,eP „r 91M 2 
"T-IIKlll^Ce) 



d-1 

E 



< |(/,<) + ( ff ,<)r B | 

2 \ 



L 2 (e 



^ ll U hllL2( rj? ) 



£ £ (t ' 



i =1 eef' 



L 2 (e 



< \(Lul) + (g,ul)r R \ 



From (|3.5p and (|3 . 6[) we can bound |it' ^ h and the jumps in terms of ||^|| i2 (f2)- 
In order to get the desired a priori estimates, we need to derive a reverse inequal- 
ity whose coefficients can be controlled. Such a reverse inequality, which is often 
difficult to get under practical mesh constraints, and stability estimates for it? will 
be derived next. 



Theorem 3.1. Let u q h £ V? solve (|2.15p and suppose (3\. e > 0,7o, e ,- ■ • , 7g.e > 0. 
Then 



(3-7) IKIIz-W + jKlliA. + KW^ + rS E W^&ie) 



i( £ c D (fc 2 K|£ a(e) + ||V4|| 2 L2(e) ))" < C«> iq M{f,g), 



where 



(3.8) M(/, fl ):=||/|| L2(n) + y| L 2 (rR) , 
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P 



1 fpk 2 h 2 +p 5 p 
- max I 1 



fc 2 eef^ V 70,e hi 



max 



h e 0<j<q-l Y 7j + l 



7j,e 



/?,„ 



+ i 



K V 7l,e 



,29+3 



// 2 



ifq<p> 



i fpk 2 h 2 e + P 5 p rj— 

— max — h — max , / 

k e ^' h \ lo, e hj h e o<j< q -i y 7 i+1 



Proof. We divide the proof into three steps. 

S^ep i: Derivation of a representation identity for ||u^||i,2(n)- Define w/j £ £ by 
f/ilif = a • Vu?|if for every if e 7^. Since is a polynomial of degree no more 
than p on K, hence, Vh €E VF. Using this Vh as a test function in (|2.15[) and taking 
the real part of the resulted equation we get 

(3.10) - k 2 Re(u q h , Vh) = Re((/, Vh) + (g, v h )r R - a p h (u q h ,v h ) - ik (u g h , v h ) rR ) ■ 
It follows from ([3TT]l . ([3~4|> . and (|3~10|) that (compare with (4.11) of [27]) 



(3.11) 2A 2 K||* 2(n) =k 2 £ I a-n K \ul\ 2 + (d-2)Re((f,ul) + (g,ul)r R 

K£T h JdK 

~ a h( u l> u h)) +2Re((/,u h ) + (g,v h )r R - a p h {u q h ,v h ) - ik (u q h ,v h ) TR ) 
= k 2 Y, I a-n K \u q h \ 2 + (d-2)Re((f,u q h ) + (g,u q h )r n ) 

KeT h JdK 

+ 2 Re((/, v h ) + (g, v h )r R ) + 2fc Im «> ^) r „ 
- Yl ((d-2)\\\7u q h \\ 2 L2iK) +2Re(\7ul\7v h )K 



KeT h 



+2 S( ( ,_ 2)Re ({g}, K1 ) /Re ({|4}, K1 

+ Re^KJ,||^|^ ^+2Im(L 1 «,^)+^J 7 « ! v h ) 



■ 1 2 I 1 2 — 

Using the identity a — 6 = Re(a + 6) (a — 6) we have 
(3.12) Y / oc-n K \u g h \ 2 =2 ^ Re(a-n e {u q h } ,[u q h ]) e + (a ■ n n ,\u 9 h \ 



2 ) 

/an • 
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Using the identity again followed by the Rellich identity (13. 2|) we get (compare with 
(4.13) of [27]) 

(3.13) £ ((d-2)||V<||^ (K) +2Re(V<,V^)K) = £ f a-n K \Vu q h \ 2 

K£T h KeT h JdK 

= 2^R e (a-n e {V<},[V<]) e + £ (a ■ n e , |V<| 2 ) e 

= 2 ^ Re(a-n e {Vul},[Vul]) e +J2 (a • n e , |V<| 2 > e 

- J] (a-n e ,|V<| 2 ) e . 

Plugging (gTTll) and (f3~TU)) into (f3~TTj) gives (compare with (4.15) of [17]) 

(3.14) 2fc 2 K||^ 2(n) 

= (d - 2) Re((/, 4) + (s, <)r R ) + 2 Re((/, « h ) + (g, v h )r R ) 
+ 2k 2 5] Re (a ■ n e {««} , [<]) e + fc 2 (a • no, |<| 2 > an 



+ 2fcIm«, Wft ) rR - 5] (a.n e ,|V<| 2 ) e + 53 <a-n e ,|V<| 2 ) 

- E He({M}, K1 ) +2 (,- 1 , E ^({gl.K, 

+ 2 53 Re^-(a-r le {V<},[V<]) e + ^{g} ) M^J 

+ 2 £ Re/[<],{|^)\ +2Im(L 1 «,^) + £j J -K,^) 

Step 2: Derivation of a reverse inequality. Our task now is to estimate each 
term on the right-hand side of (|3.14p . Since the terms on the first four lines can 
be bounded in the exactly same way as done in [27j . we omit their derivations and 
only give the final results here for the reader's convenience. 

(3.15) 2Re((f,v h ) + (g,v h )r R ) 

<CM{f,gf + \\ul\l h+ C -^Yl 



4 II v "Tilli^e) • 



(3.16) 2fc 2 53Re(a-n e K},KJ) c 



^ fc2 II ? II 2 i n P fc2 7o,eP I, r 91 ,,2 

-2 |L,9|| 2 



(3.17) fc 2 (a • «„, |<| 2 ) ao < Ck 2 K|| 2 2(rj?) + £ fc 2 (a • n e , | u «| 2 > { 
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(3.18) 2klm(ul,v h ) rR - E (a-n e ,\Vu q h \ 2 ) e 



ees! 



(3.19) 2(d- 1) E Re 



l <9n e J 



<l\<\l h +c E 



P 70,eP ||r gi|| 2 
— ^HWII^e)- 



The extra work here is to estimate the terms on the last two lines in (|3 . 14[) . For 
an edge/face e G El , let f2 e denote the set of element (s) in Th containing e as one 
edge/face. By (|3.3p we obtain 

(3.20) 2 £ Re f- (a • n e {V<} , [V<]) e + 



\ <9n e J 



d-l 

=2 E E Re 

eeS! D 1=1 



a ■ t„ 



dn e 



a ■ n, 



1 dri J 



< E pV E llv<ll i2W IIK]|| i3(e) 



ees! 



Ken c 



<Ik\Ih+ c E 



P 7b,eP „r g1 „2 



The first term on the line six of (|3.14[) is bounded as follows (compare with (4.20) 
of [23): 



(3.2i) 2 E Re ([<],{Sr|) £ E ^iK]|| L2(e) E 

eef^ X 1 e e e£, JD Aesi c 



L2(if) 



< E ^ 



4UHz, 2 (e) E H^ n 0z 2 (lT) 

A'eo e 



<~ 8 K\i, h +c E iikjw 



The penalty term £i(-, •) is estimated as follows. Recall that Vh\K — ® ■ ^ u 1\k 
with a — x — xq 1 for each K £ 7j t . Noting that 



(3 22) ** = *£+a-v(^ 



9rf 



+ a ■ n ( 



d-l 



dri V dn. 



) + E«- 



d rdu q h \ 



KKd-l, 
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by the definition of L\(-, ■) and Lemma l3~2l we get 



d-l 



(3.23) 21mLi(u q h ,v h ) =2Im E E 



Pl,eP 



a ■ n e 



d / du q h n 



d-i 



< 



< 



2i- EE 

eee° t=i 

d-l 

+ ^EE 

d-l 

+ ^EE 

d-l 

21- EE 

edSP 1=1 



d-l 

+ E 

m=l 

Pl.eP 

h e 

Pl,eP 3 

Pl,eP 3 
hi 



du q h 



9uj d fdu q h 
dri ' 3t{ V dn e 



9ri 



dri 



L 2 (e) 



du q h 



LHe) 



Pl,eP 



au H h a / au H h 



d-l 



P 3 (3l,e / 7l,e/le 



<9n e 



9ri 



L2(e), 



+ ^EE 



P 2 Pl,eP 
he K 



'dul' 

dri 



L 2 (e 



We remark that Im L\(u\,Vh) = when p = q = 1. 

Next we estimate the penalty terms Jjiu^Vh)- Since u q h and vh are piecewise 
polynomials of degree p in general and those terms contain jumps of high order 
normal derivatives, it is quite delicate to control those terms as shown below. 

By direct calculations we get that on each edge/face e of K 6 7/, 



(3.24) ^=j-%+a^(^i) 

d*ui d j+1 u q h m 9 (d j u q \ 
= 3 1+a-rie -pr+> oi-tV 1 - — ( M, l<7<p-l, 

/o«^ ^"h 5P < ^fd p u q \ d p u q h 



Here we have used the fact that (p + l)th order derivatives of u q h is zero because u q h 
is a polynomial of degree at most p. 



14 



XIAOBING FENG AND HAIJUN WU 



For j = 1, 2, • • • q - 1, by (f3T24j) we have 



(3.26) 21m Jjiulvh) = 21m ^ 7j - e f 



2j-l 



dh 



dnl 



dnl 



+i 



E Q ' 



dnl 



( dj <)]\ 
\ dnl )\/ t 



dVi 



S E 7i,. 



2J-1 



L=(e) 



-1^9 



dnl +1 



< JL / 7?> 



L2(e) 

2i-i 



9ni 



L2(e) 



23+1 



(9rie 



+7j+l,e 



2j-l 



9^ 



£ 2 (e), 

2 



L=(e) 



If g < p, then, from Lemma 13.21 and the inequality 
we have 



&><p 



dnl 



L 2 (dK) 



m{Ky 



(3.27) Im2J 9 (u« > « A )< £ 7 „ 



2g-l 



dni 



* E 7,.. ( 7 



9+2 



£=>(e) 

0«< 



Vh 



m(K) 



dnl 



i2 ( e ) *ren e 



ft-? 





9 9 < 




dnl 



L 2 (e) 



If q = p, (I3.25| and the definition of J p (-, •) immediately imply that (compare with 
(4.14) of [27]) 

(3.28) 21m J P (u q h ,v h ) = 2ImJ p (u 9 h ,u q h ) = 0. 

The estimate for ImJo(«',Uh) is similar to (13.26|) . so we get 



(3.29) 2ImJo«,^) = 2Ini £ ^ ( KJ 



hi > 



ee£. 



'dn e + ^ a T *dri 
J= i 
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lO.eP 



< 



2Im E 



ee£f 



dn e / 



+^E 



70,e£> 



l hl\\L 2 (e) 



du q h 



lO.eP ... g 
—. U 



dn e 



< 



2Im E 



10,eP 



a ■ n e u\ 



q 9u h \ 
' dn e / 



^E 



ee£! 



P lO.eP „r g 1M 2 

r-f- \\[u h }\\ LHe) 



+ ^E 



P /70,e / 70,eP 



ee£ 



\ n e V 7i, e 



i ?iJlli 2 (e) 



We also need the following estimate (compare with (4.22) of [27] ) 
(3.30) £ (V(a.n e ,K| 2 ) e + (a.n e ,|V<| 2 ) e 



L 2 (e 



d-1 



£=1 



2^ h . Lm \ a n - d T l'dTi\dn e 
dn e I 



Im ( a ■ n e u\, 



eest 



<-J2 (a-^,fc 2 |<| 2 + |V<| 2 -2^ 



01, eP 









dri 


dri 


\dn e J 



- 2 



70,eP, 9 



c 2^ —f^l^-f— 



dul 



dri 



L 2 (e 



2 

i 2 (e) 



where we have used the inverse inequality and the assumption that D is star-shaped 
to derive the last inequality. 

Step 3: Finishing up. We only prove the case of q = p since the proof for q < p 
is the same except using (pHTf instead of (pT2"g|) . Substituting (pHS)) 
(with q = p) and (|3~TS)) - (l3~Tg|) into (l3~T4")l . and using (pPO"]) we obtain 



2fc 2 |k p H 2 



7illL 2 (f2) 



<(d-2)Re((/,<) + (<?,<>rj + ^Af(/,s) 2 + T Kll£ a(n) 



l h\l,h 
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CD 



■2|L,P||2, 



*felli a (e) "T" j ll^ n h.lli 2 (e) 



P\\2 



+ ^E 



P(fc 2 + 1) P 5 70,eP 

70, e 70,eAe ^ 



C J2 



70, eP , P 



P \ 10,eP 



lQ,eh 2 e h e 



l hi\\L 2 (e) 



d-1 



ee££ 



du p h 



d-l 



F 5 //3l,e / 7l,e^e 



^ E E 



P-1 



P 2 ffl.eP 
ft fi /ifi 



p-1 

i=o eeef e V 7j+ 



P 



L 2 (e) 



L 2 (e) 



L 2 (e 



3»i 



2j-l 



dnl 



+7j+l,e [ — 



2j+l 



di+ l u p h 



dnl 



L 2 (e), 



Therefore, it follows from Lemma and (13.91) that 



2/ fc2|U,P||2 



CD 



■2|U,P||2 



*fcllL 2 (e) 2 



Jliv<ni 2(e) 



< CM(f,g) 2 + — \\u P X 2{Q) + Ck 2 C sta , p + ( 3 ,<)r„| 
<Cfe 2 C B 2 ta , p M(/ !5 ) 2 + ^-K|| 2 2(n) , 



where we have used the following inequality, which is a consequence of (|3.6j) . 



k 2 \\ui\\\ 2(rR) <e\\uir L2m +M(f,gf 



-F/P-DISCONTINUOUS GALERKIN METHODS FOR THE HELMHOLTZ EQUATION 17 



to derive the last inequality. M(f,g) and C s ta, P are defined by p.8|) and (13. 9|) 
respectively. Hence, 



?llr a , m +^|tt?L .+T- CO, V IIV^I 



x l p 1 









l"fcllL»(n) + l\ U h\l,h + k\ Cn i 2^ H V %Hi=(e) 



1 

+ l( E ^( fc2 KH£'(e) + HV<ll! 2(e) )) 2 < C sta , p M(f,g), 

which together with (|3.6|) gives (|3.7|) . The proof is completed. □ 

As (|2.15|) can be written as a linear system, an immediate consequence of the 
above stability estimate is the following well-posedness theorem for (|2.15|) . 

Theorem 3.2. For k > and h > 0, the hp-IPDG method (|2T5|) has a unique 
solution u q h provided that 7o,e>7i ei ' ' '7g,e > and /3i, e > 0. 

Next we consider the case of quasi-uniform meshes. Note that large penalty 
parameters ("fj,j > 1) for jumps of normal derivatives may cause a large interpo- 
lation error in the norm h , and hence, may pollute the error estimates of the 
IPDG solution (see Section HJ. It is interested to minimize the stability constant 
Csta, 9 under the constraints of /?i >e > and + $Z' =1 P 2j_1 7.?> ~ !• We have the 
following consequence of Thcorcm l3.ll The proof is straightforward and is omitted. 



Theorem 3.3. Let h = max/i e . Suppose the mesh Th is quasi-uniform, that is 
h e ~ h. Suppose k > 1 and kh < 1. Assume that 7^ ~ 7j,j = 0, 1, ••• ,q, 

E|=iP 2j ' _1 7j < 1, < /?i, e < ^7o, <md tkrf 



J 7o~_p3ft, 3, 7j>p 3/i3 7i _ x ^ ifq = P, 

[ 7o ~ min jpT£"/riTT,pi/i-t j , 7j > (^) 7i-i> 7g ^ 7oP ^ l f 1 < P- 
Suppose D — if q > 2 or q = p = 1 . Tften 

IKIli^) + £ IKIIl,h,9 ~ C *t a ,qM{f,g), 

where 

p3 



Csta.o 



fc 2 /l3 



pap 9+1 



,,,41 „ g + 2 



ifq = p, 

if q< p. 



It is clear that, in the above theorem, jo — P 3 h 3 if q = p or q > 2, 70 — 

3g+l q_ 

p 9+1 ft 9+1 otherwise. We conclude this section by several remarks. 

Remark 3.2. (a) The hp-IPDG method (|2.15|) is well-posed for all h, k > pro- 
vided that all penalty parameters are positive. As a comparison, we recall that the 
standard finite element method is well-posed only if mesh size h satisfies a constraint 
h = 0{k~ p ) for some p > 1, hence, the existence is only guaranteed for very small 
mesh size h when wave number k is large. 

(b). It is well known that [3 HI El EH SQl [45] symmetric IPDG methods for 
coercive elliptic and parabolic PDEs often require the penalty parameter 70, e is suf- 
ficiently large to guarantee the well-posedness of numerical solutions, and the low 
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bound for 7o, e is theoretically hard to determine and is also problem- dependent. 
However, this is no issue for scheme (12.15P , which solves the ( indefinite ) Helmholtz 
equation, because they are well-posed for all 7o, e > 0. 

(c) In the linear element case of q = p — 1, a better estimate is obtained in |27] 
with the help of the penalty term L±. Unfortunately the penalty term L\ does not 
help too much for higher order elements. 

(d) . The stability estimates will be improved greatly when k^h^p -1 < Co, where 
Co is some constant independent of k, h, p, and the penalty parameters (see Theo- 
rem \5.1\) . The above estimates are only interested when kh < 1 and k z h 2 p^ 1 > 1. 

4. Error estimates 

In this section, we derive the error estimates of the solutions of scheme (|2.15[) . 
This will be done in two steps. First, we introduce elliptic projections of the PDE 
solution u and derive error estimates for the projections. We note that such a 
result also has an independent interest. Second, we bound the error between the 
projections and the IPDG solutions by making use of the stability results obtained 
in Section [3l In this and the next section, we assume that the mesh 7^ is quasi- 
uniform, that 7j !e ~ 7j > for j = 0, 1, ••• ,q, and that [3i >e ~ f}\ > 0. Let 
h = max/i e . For simplicity, we also assume that the mesh 7/j is conforming, that is, 
Th contains no hanging nodes, since the parallel results for non-conforming meshes 
can be derived in a similar way. 

In this and the next section the following assumption (cf. Theorem I2.1|) is re- 
quired for some results. Suppose problem ()l.ip - (|1.3|) is IP-regular and 

(4-1) hWns^^k^M^g), 

where s > 2 is an integer and 

M s (f,g) = ||/|| H „,ax {0 , s -2} (n) + ||5lU 2 (r„)- 

Note that the assumption is proved to hold for s = 2 and Ma(/, g) = M(f,g). 

4.1. Elliptic projection and its error estimates. For any w G E q fl Hp D (£l) n 
H? oc (Q), we define its elliptic projection € V£ by 

(4.2) al(wl,v h ) + ik(w q h ,v h ) rR =a q h {w,v h )+ik(w,v h ) ra ~iv h G V% . 

In other words, is an IPDG approximation to the solution w of the following 
(complex-valued) Poisson problem: 



dw 
dn r , 



-Aw 


= F 


in f2, 


+ ikw 


= i> 


on Tfi 


w 


= 


on r D 



for some given functions F and ip which are determined by ui. 

Before estimating the projection error, we state the following continuity and 
coercivity properties for the sesquilinear form a q h {-, ■). Since they follow easily from 
(|2.4p ~ (|2.13p , so we omit their proofs to save space. 

Lemma 4.1. For any v £ E q and w G E q n Hy D> the mesh- dependent sesquilinear 
form a^(-, •) satisfies 

(4-3) l<KHI,KK«)l<lkll Vl ,Jlkllli,M. 
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In addition, for any < e < 1, there exists a positive constant c £ independent of k, 
h, p, and the penalty parameters such that 

(4.4) Real(v h ,v h )+(l-e + ^)lma q h (v h ,v h )>(l-e)\\v h \\ 2 lhg Vv h eV%. 

\ 70 J • ' H 

Let u be the solution of problem 11 . 1|) — ()1 .3|) and be its elliptic projection 
denned as above. Define r\ :— u — u q h . Then (14. 2\ immediately implies the following 
Galerkin orthogonality: 

(4.5) a q h (r,,v h )+ik(ri,v h ) rR =0 V^eV^f. 

Lemma 4.2. Suppose problem (|l.l[) - (|1.3p is fj niax {i+ 1 - 2 } -regular. Then there hold 
the following estimates: 



(4.6) ||u - u q h \\ l h q + VXk \\u - u q h \\ L2{TR) 



< inf 

z h ev*nHl in) 



(\\\\U - Zh\\\l,h,q + VAfc ||u - z h \\ L 2 (rii ^j 



(4.7) || u -^|| i2(n) </ l (l + -L + ^+M 



loP P 

x M i (x\\\u-z h \\\ 1Jhq + VXk\\u~z h \\ L2{rR) ) 



where A := 1 + — and 71 = if a = 0. 



Proof. For any z/j S H -ffp D (f2), let rjh = u q h — Zh- From rjh + n ~ u — z% and 
(14.51). we have 



(4.8) aliVh^h) + ik {VhiVh)r R = a K u ~ z >» Vh) + ik(u- z h ,r) h ) TR ■ 

Take e = \ in (|4.4p and assume without loss of generality that ci > i- It follows 
from (I4~4l) and (l4~8l) that 



1 2 /l C -P\ 

2 ll%lli,h,g < Re o£ + (_ + -2_J tena q h (r)h,rjh) 

/ 1 cip\ 

= Re(a 9 h (u- z h ,T) h ) + ik{u- z h ,n h ) rR ) - ^- + -^-)k(rj h ,r) h ) rR 
+ (- + Im (a£(zt - 2ft, r? h ) + ifc (u - z h ,i] h ) rR ) 

< CX {\\Vh\\ hhtq \\\u - Z h \\\l,h,q +k\\u- Zhf L *F R )) - -J- \\Vh\\ 2 L 2 {rR) • 

Therefore, 

(4-9) IMIi.m + Afc ll%lli2(r B ) < a2 III u - + Afc ll u " z hlli,=(r R ) 



which together with r\ = u — zu — Tjh yields 

To show (|4.7p . we use the Nitsche's duality argument (cf. [HHH]). Consider the 
following auxiliary problem: 

(4.10) -Aw = 77 in (1, 

l/sw = on T R , 

on 



w = on r 
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It can be shown that w satisfies 

( 4 - n ) Hff>(n) ~ Nli»(n) • 

Let w\ be the continuous linear finite element interpolant of w on 7^. From (|4.5p . 

aliv, $1) + ik {v, ™h) rR = a l (V, + ifc {v, ™1) Tr = °- 
Testing the conjugated (|4.10|) by rj and using the above orthogonality we get 

IMIi=>(n) = a Uv,w) + ik (v,w) Tr 

= aUv, w - w\) + ifc (77, w - w\) Vr 

^ IMIi,h,i Hl w ~ ^IHw + k WvWl'Fr) \\ w ~ ™fclL=(r fl ) • 

£ IMIi,M ( x + ^ + j) 2 HH^o) + k Mi*?*) hi Hn*(ny 
which together with (|4.6|) and (|4. gives (|4.7p . The proof is completed. □ 
We have the following lemma that gives approximation properties of the space 

v^nif^(O). 

Lemma 4.3. 

(i) Let /i = min {p + l,s} and q < \x. Suppose u G H S (Q) n ifp (f2). Then 
there exists G V^ 1 PI i/p D (51) smc/i i/iat 

ft/*- 3 

(4.12) ||u - %||i2 ( r R ) < -^r IMIff'(n) . 

p 2 

(4.13) |||« - fi fc ||| 1>M < (l + JL + J2p 2 'S) 2 ~~~T IMI^(n) ■ 

7o j=1 P 

(ii) Suppose u£ if max {« +1 > 2 > (O^i/j^ (ft). Then there exists u h G V£nH£ D (Q) 
such that l|4.12p holds with s = 2 and 



9 It 

(4.14) |||u - u A |||i,j,,, < (l + — +P7i + J2 p2j ^^ 



70 — / P 



where 71 = if q = 0. 



Proof. The following Zip approximation properties are well-known for the hp finite 
element functions (cf. [51 IBTil I3"T] ): 

• There exists Uh G V? such that, for j = 0, 1, • • • , s, 



"II 

pa 



(4.15) \\u-u h \\ HH r h ) ■= ( H IN-^llL'CK)) 2 ~^JH"ii// ,<>> 

• There exists Uh G n Hp such that 

h^~ j 

(4-16) II^-%IIh3(o) ^ — " IMIffa(n) > J = 0,1. 
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Here the invisible constants in the above two inequalities depending on s but in- 
dependent of k, h, p, and the penalty parameters. Then (I4.12| follows from (|4.16p 
and the trace inequality. 

It follows from the inverse inequalities in Lemma 13.21 that . for 1 < j < q + 1, 



(4.17) 



\u - u h \ 



HHT h 



< h - u h \\ HHTh) + \\u h - Uh\\ H3{Th ) 



< 



u - u h \ 



p 



Hi(T h ) + frj-i W U h - u h\\ H 1 (n) 



~ p s-2j + l \\ U Wh'((1) 



Therefore, by the following local trace inequality 



we have 
(4.18) 



E E 

KeT h eCdK 

s.h- 1 



d 3 {u - u h ) 



L3(e) 



u ^h\\Hi(T h ) + \\ u ^' l llffJ(T h ) ll M "fc|lffi+ 1 (75 l ) 



< 



h lH-2j-l 



P 



2s-ij 



Noting that Uh is continuous, we have from (|4.16[) and (|4. 1 8[) . 



(4.19) |||u - u h \\\l, h>q =\u- u h \\ h + Tj,e [j 



3=1 e££f 



d^u - uh) 



E 



<9(u - u h ) 



L 2 (e) 



< 



P 2," 1 „ „2 



(O) 



dni 



me) 



That is, (|4~T3]) holds. 

(I4.14| can be proved similarly as above. It is clear that (|4.16p and (|4.17p hold 
with s = 2 and (14.18P holds with s = 2 and j — 1, that is, 



(4.20) 
(4.21) 

(4.22) 



Iff 2 (a) 



|u-u/i||jr»(T h ) ^Pll M ll^ (n) , 



E E 

KGT h eCdK 



d(u - U h ) 



dn e 



<h\ 



j = 0, 1, 



\H 2 (Q.) 
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Since u e H q+1 (Q), we have from Lemma I3~2l and (I4.21|) that, for 2 < j < q, 



(4.23) 



E 

eeel 



d J (u - u h ) 



dni 



< 



< 



E E 

KeT h eCdK 
,2 s 2j-3 



L 2 (e 

d 3 Uh 



E 

eesl 



d 3 iih 
dni 



L 2 (e 



di 



< 



\u h 



H (Th) 



LHe) 

< 



P || - ,|2 

-T \\ u h\\ H j(T h ) 



fP\ 2j - 3 2j-2 II ||2 



(Th) 



Now dHHl) follows by combining the equality in flU?]) and P~2"0]) - (|3~2"3")) . This 
completes the proof of the lemma. □ 



By combining Lemma 
the projection error. 



and Lemma I4.3f i) we have the following estimates for 



Lemma 4.4. Let /i = min{p+ l,s} and q < u. Suppose problem - (jl.3[) is 
H s -regular and (14. ip holds. Then there hold the following estimates: 



(4.24) 

(4.25) 

where 



"fclll.h,. 



f- VAfc||n-ti^|| i2(rjj) < C c 
<f C^M.t/.j)— , 



; q M s {f,g) 



p s-l 



Ccrr,q -~ 



a 



crr,q 



A(l + ^ 

7o 

1 

7o P 



E^~S 



:= 1 



— + kh 



a 



crr,g ; 



A:=l + ^. 

7o 



Remark 4.1. TTie requirement q < s in the lemma is clear since the projection u h 
is not defined for q > s. However, for q < s, \\u — u q h \\ 1 h can be bounded without 
using full regularity of u, and such a bound is also useful (see Lemma \5.1\) . 



4.2. Error estimates for 



In this subsection we shall derive error estimates 



for scheme (I2.15[) . This will be done by exploiting the linearity of the Hclmholtz 
equation and making use of the stability estimates derived in Theorem 13.11 and the 
projection error estimates established in Lemma 14.41 

Let u and u q h denote the solutions of (|l.l|) - l|1.3p and (|2.15|) . respectively. Assume 
that u S H s (n) with s > q + 1, then f2~T4")) holds for v h E V%. Dehne the error 



function th 
equation: 

(4.26) 



Subtracting (|2.15p from (|2.14p yields the following error 



a q h (e h ,v h ) - k 2 (e h ,v h ) + ik(e h ,v h )r R = \fv h G V%. 
Let ui be the elliptic projection of u as defined in the previous subsection. Write 



eh = rj — £ with n 



From (|4~2^| and (113)) we get 



(4.27) al(£,v h ) ~ k 2 (£,v h ) +ik{£,v h )r R = a q h (T],v h ) - k 2 (n,v h ) + ik(r},v h ) 

= -k 2 (r],v h ) Vv h eV?. 
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The above equation implies that £ £ V£ is the solution of scheme (|2.15| with sources 
terms / = —k 2 r\ and g = 0. Then an application of Theorem 13.11 and Lemma 14^41 
immediately gives 

Lemma 4.5. £ = u q h — u q h satisfies the following estimate: 

1 ^ h s h IJ - 

(4.28) ||C|| La(n) + ^ ||£|| lih ,, < C^, q kpC m , q M a (f,g) — . 

We are ready to state our error estimate results for scheme (|2.15|) . which follows 
from Lemma I4.5[ Lemma 14.41 and an application of the triangle inequality. 

Theorem 4.1. Lef u and u q h denote the solutions of (|l.l|l - (jl.3|) and (12.15|) . re- 

spectively. Let fi = min{p + l,s} and q < fi. Assume assumption (14. II) holds. 
Then 

/ ~ \ k s ~ 1 h fJ — 1 

(4.29) ||u - u«|| lifcifl < (C err , g + C sta , g fc 3 /iC ori , g ) M s (/,<?) — , 

- /I \ k s h>* 

(4.30) ||« - u«|| ia(n) < P C m , q (- + C sta , 9 /c) M s (f,g) — -. 

Remark 4.2. q < s is required in the theorem because \\u — u^]^ h is not defined 
for q > s. However, we note that the hp-IPDG solution u q h is always well-defined 
regardless the regularity of underlying PDE solution u. For q < s, \\u — u q i \\ 1 h 
can be bounded without using full regularity of u, and such a bound is also useful 
(see Lemma\5., 



By combining Theorem 14 . 1 1 and Theorem 13 . 31 we have the following theorem that 
gives the best convergence order so far that we can obtain theoretically for the 
method (|2.15[) under the mesh condition fc 3 /i 2 p _1 > 1 (cf. Theorem 15. ip . 



Theorem 4.2. Under the assumptions of Theorem \3.S\ and \4.1\ we have 
(4.31) \\u-ul\\ lhq + k\\u-ul\\ LHn) 

< f kpih-i M s (f,g) fc'^r 1 ifq=P, 
~ \ fc max |pf h-^,p^hsh\ M s (f,g) ^"ffi -1 ifq<P- 

Proof. The proof is obvious since C crr ,q — C cr r,<j — 1. Q 

Remark 4.3. (a) Estimates ()4.29|) (j4.31|l are so-called preasymptotic error esti- 
mates which are suboptimal in h and k. They can be improved to optimal order 
when k^h 2 p~ x < Co, where Co is some constant independent of k, h, p, and the 
penalty parameters (see Theorem \5.1\) . The second term on the right-hand side of 
(I4.29|) is called a pollution term for \\u — u q h \\ 1 h q . 

(b) Theorem \4-£\ shows that \\u — u^W^ h q + k ||w — u^Wl 2 ^) ^ */ 1 = P an d 
p^~ s k s h^~^ — > 0, or if q = 1 < p and p i ~ s k s h ll ~i — > 0. 

5. Stability-error iterative improvement 

In this section we derive some improved optimal order stability and error esti- 
mates for the hp-IPDG solution under the mesh condition that k 3 h 2 p~ x < Co by 
using a stability-error iterative procedure, where Co is some constant independent 
of k, h, p, and the penalty parameters (see Theorem 15. ip . 
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By combining Lemma 14.21 and Lemma 14.3( h) , we have the following estimates 
for the projection error when only the £f 2 -norm of the solution u is allowed in the 
error bound. 

Lemma 5.1. Suppose problem (jl.l|l - (|1.3j) is H max { q+1 ' 2 } -regular. Then there hold 
the following estimates: 

k h 



(5.1) 


ll u - fi hlli,fc, g ' 


(5.2) 


\\ u -KWv(a) 


where 






Corr,2,(j '•= A(l 




C er r,2,g '•= ( 1 ^ 



P ~ , . k h 



2 h 2 



~ k ° IT '" q p 2 ' 

P kh 
7o Ap 

— + ^ + khY C m , 2 , q , \:=l + 2-. 
70 p p / 70 

By a similar argument to that used to prove Theorem l4.H we have the following 
error bounds which only involves M(f,g). 

Lemma 5.2. Let u and u q h denote the solutions of (| 1 . 1 [) (| 1 .3[) and (]2 . 1 5[) . respec- 
tively. Suppose u G H max {« +1 ' 2 >(n) n i?r D (fi). T/ieri 

(5.3) < (Corr,2,,+C B t a ,,fc 3 /lC orr ,2,,)M(/, fl )— , 



fc 2 /* 2 

2~ • 



(5.4) ||«-<||^ (n) <Pa»,2 l9 (^ + C llte , 9 fc)M(/, fl )^ 

We are now ready to state our final main theorem of this paper. 

Theorem 5.1. Let u and u q h denote the solutions of (|l.l|) - (|1.3j) and (I2.15|) . re- 
spectively. Suppose u e H max ^ +1 ^(fl) n H£ D (il). Assume that k>l,kh<\, 
and that P7o _1 + 2?=i P ij 1^ 1- 27ien i/iere exists a constant Cq > 0, which 
is independent of k, h, p, and the penalty parameters, such that if k 3 'h 2 p~ 1 < Cq, 
then the following stability estimates hold: 

(5-5) K\\i, h , q ZM(f,g), 
(5-6) Kll^n) 5 ). 



Moreover, under the assumption (|4.ip . f/iere /ioW i/ie following error estimates: 
(5-7) ||u - < (1 + k 2 h) M s (f, g) -± , 

(5-8) \\u~ul\\ L2m < P M s (f,g)^. 

Proof. We only prove (|5.5[) since (|5.6[) can be proved similarly and (|5.T[) — (|5.8p fol- 
low from the improved stability estimates and the argument used in the proof of 
Theorem 14.11 

From Theorem 13. II we have 

(5-9) \\ul\\ lhiq <kC Bta , q M(f,g), 
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where C s ta,g is defined in (13.91) . From Lemma I5T21 we have 

(5.10) h-u q h \\ hhq < (c e „, 2tq +kC statq k 2 hd e „, 2 , q ) — M(f,g). 

Now it follows from Theorem 12. II and the triangle inequality that 

( 5 - n ) IKIIi.m - W u hh, q + W U - u l\\i,h, q = \ u \hh + h- K\\i, h , q 

( kh k 3 h 2 ^ \ 

< 1 + C* 0IT .2,g — + C err ,2,g k C sta .g Af (/, g). 

V p p J 

Repeating the above process yields that there exists a constant C\ independent of 
k, h, p, and the penalty parameters, and a sequence of positive numbers Aj such 
that 

(5-12) H\\i, h , q <*jM(f,g), 
with 

kh -» k 3 h 2 

A ~ kC sta q , Aj = Ci(l + C crr 2,o — ) + Ci C crr 2,0 j ' — 1;2, • • • . 

p p 

A simple calculation yields that if C\ Con, 2,0 k 3 h 2 < Op for some positive constant 
< 1 then 

,. . Cl(p + Ccrr, 2,o kh) 

hm Aj = — — — , 

P~Cx C crr ,2,o k 3 h 2 

which implies (|5.5p by noting that C eTI> 2 >q , C er r,2,g ^ 1 and that C cir .2, q kh < 

Note that the stability estimates in (|5.5p and (|5.6p are of the same order as the 
PDE stability estimates given in Theorem 12. II 
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